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Abstract 

This paper is concerned with the inverse problem of scattering of time-harmonic acoustic 
waves by an inhomogeneous penetrable obstacle in a piecewise homogeneous medium. The 
well-posedness of the direct problem is first established by using the integral equation method. 
We then proceed to establish two tools that play an important role for the inverse problem: 
one is a mixed reciprocity relation and the other is a priori estimates of the solution on 
some part of the interfaces between the layered media. For the inverse problem, we prove in 
this paper that both the penetrable interfaces and the possible inside inhomogeneity can be 
uniquely determined from a knowledge of the far field pattern for incident plane waves. 

Keywords: Uniqueness, piecewise homogeneous medium, penetrable obstacle, unique 
continuation principle, Holmgren's uniqueness theorem, inverse scattering. 

1 Introduction 

In this paper, we consider the scattering of time-harmonic acoustic plane waves by an inhomoge- 
neous, penetrable obstacle in a piecewise homogeneous surrounding medium. In many practical 
applications, the background medium might not be homogeneous and then may be modeled as 
a layered medium. We might think of a problem in medical imaging where we have a damaged 
tissue (which can be modeled as a penetrable obstacle) inside the human body. It is clear that 
the human body is not a homogeneous structure and may be modeled as a piecewise homoge- 
neous medium. Therefore, one possible model would then be a penetrable obstacle buried in a 
piecewise homogeneous medium. 

For simplicity, and without loss of generality, in this paper we restrict ourself to the case 
where the obstacle is buried in a two-layered piecewise homogeneous medium, as shown in Figure 
[TJ Precisely, let C be an open bounded region with a boundary Si such that the 
background M^\f22 is divided by means of a closed surface Sq into two connected domains JIq 
and Oi. Here, f^o is a unbounded homogeneous medium, J7i is a bounded homogeneous medium 
and Q.2 is a penetrable obstacle. Let denote the complement of VLq, that is, Q. := M.^\Qq. 
Choose a large ball Br centered at the origin such that C Bji and let fi/j := Bji\Q. 

The problem of scattering of time-harmonic acoustic waves in a two-layered background 
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where v is the unit outward normal to the interface Sj (j = 0, 1) and n € C^''^{0,2), < a < 1, 
is the refractive index of an inhomogeneous medium with 3?(n) > and 9(n) > 0. Here, the 
total field u = is given as the sum of the unknown scattered wave which is required to 

satisfy the Sommerfeld radiation condition (II. 6p and the incident plane wave = e*'^"^ '^, kj is 
the positive wave number given by kj = ujj/cj in terms of the frequency loj and the sound speed 
Cj in the corresponding medium (j = 0,1,2). The distinct wave numbers kj (j = 0,1,2) 
correspond to the fact that the medium consists of several physically different materials. On the 
interfaces Sq and Si, the so-called "transmission conditions" (jl.4p and (II. Sh with two positive 
constants Aq and Ai are imposed, respectively, which represent the continuity of the medium 
and equilibrium of the forces acting on them. 




Figure 1: Scattering in a two-layered background medium 



The direct problem is to look for a set of functions {u,v,w) satisfying (jl.ip - (jl.6p . We will 
establish the well-posedness of the direct problem, employing the integral equation method in 
the next section. We will also establish a mixed reciprocity relation and prove a priori estimates 
of the solution on some part of the interfaces Sj {j = 0, 1). These results will play an important 
role in the proof of the uniqueness results in the inverse problem. 
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It is well known that u^(x) has the following asymptotic representation 
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uniformly for all directions x := where the function d) defined on the unit sphere S 

is known as the far field pattern with x and d denoting, respectively, the observation direction 
and the incident direction. 

The inverse problem we consider in this paper is, given the wave numbers kj (j = 0,1,2), 
the positive constants Aq, Ai and the far field pattern u°°(x, d) for all incident plane waves with 
incident direction d € S", to determine the interfaces Sj (j = 0, 1) and the refractive index n. 
Precisely, we will study the uniqueness issue and prove that the interfaces Sj {j = 0, 1) and 
the inhomogeneity n can be uniquely identified from a knowledge of the far field pattern. It 
should be remarked that the uniqueness issue for inverse problems is of theoretical interest and 
is required in order to proceed to efficient numerical methods of solutions. 

For the unique determination of an inhomogeneity with compact support in a homogeneous 
background medium, we refer the reader to Hahner [5j, Nachman [13], Novikov [13] and Ramm 
[15\ [T6j . We also refer the reader to the monographs of Colton and Kress [3] and Kirsch [7\ 
and the habilitationsschrift of Hahner [6] for a comprehensive discussion. Kirsch and Paivarinta 
[9] gave the first uniqueness result for determining the interior inhomogeneity in the case of a 
known inhomogeneous background medium. 

However, to the authors' knowledge, there are few uniqueness results for the inverse obstacle 
scattering in a piecewise homogeneous medium. In the case when the obstacle ^2 is impenetrable, 
based on a generalization of the mixed reciprocity relation, Liu, Zhang and Hu [12j proved that if 
the interface 5*0 is known then the obstacle and its physical property can be uniquely determined 
by the far field pattern. Later, motivated by Kirsch and Kress [8j and Kirsch and Paivarinta 
[9] for the uniqueness proof of determining a penetrable boundary, Liu and Zhang [11] extended 
this result to the case when the interface 5*0 is unknown and proved that the interface Sq can 
also be uniquely recovered. Note that all the results in [12^ [TT] are also available for the case 
of a multilayered medium and can be proved similarly. Thus, Liu and Zhang [11] have in fact 
proved a uniqueness result for the case when the obstacle ^2 is penetrable with a homogeneous 
interior. This result has also been proved by Athanasiadis, Ramm and Stratis [Ij and Yan [19] 
using a different method. 

In this paper, we consider the case where the obstacle is penetrable with an inhomogeneous 
interior and use the ideas in [HJ [9l [11] to prove the unique determination of the interfaces Sj 
(j = 0, 1) in Section [3l In Section [31 we will show that the refractive index n is also uniquely 
determined. To do this, we first establish a completeness result, that is, the normal derivatives 
{dw{-; d) /du, d £ S}, corresponding to incident plane waves with all directions d G 5, are 
complete in L'^{Si). Based on this result, we then establish an orthogonality relation that 
enables us to prove the unique determination of the inhomogeneity n with help of the ideas from 
the fundamental work of Sylvester and Uhlmann [17J. 

2 The direct scattering problem 

In this section we will establish the well-posedness of the direct problem via the integral equation 
method. We also prove a mixed reciprocity relation and a priori estimates on some part of the 
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interfaces Sq, Si of the solution with the help of its explicit representation in a combination of 
layer and volume potentials. We assume hereafter that k^, ki, k2, Aq, Ai are given positive 
numbers and that fc^ is not a Neumann eigenvalue of Aw + felraw = in 1^2- In this paper, we 
shall use C to denote generic constants whose values may change in different inequalities but 
always bounded away from infinity. 

Remark 2.1. The assumption on fe^ holds provided the refractive index n satisfies one of the 
following two conditions: 

(1) 9(n) > on some non-empty open subset 0* of i^2', 

(2) 1 — n is compactly supported in 0,2- 

Both cases can be proved by the unique continuation principle. 

As incident fields n*, plane wave e*^"^*^ and point source ^j{-,Zj), (zj G Qj, j = 0,1), (cf. 
(j2.8p below) are of special interest. Denote by u^{-,d) the scattered field for an incident plane 
wave u^{-,d) with incident direction d € S and by u°°{-,d) the corresponding far field pattern. 
The scattered field for an incident point source zj) with source point Zj G Qj is denoted by 
Zj) and the corresponding far field pattern by <I>°°(-, Zj) (j = 0, 1). 

The direct problem is to look for a set of functions u G C^(Oo) ^ C^'"(Qo)) ^ C^(f^i) H 
Ci'"(ni) and w e c^{n2) n Ci'"(n2) 
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p G C"'^'"(S'i) and q G C^'"{Si) are given functions 


from 



where / G C^''^{Sq), g G C°'°(^ 
Holder spaces with exponent < a < 1. 

Remark 2.2. For the scattering problem, if the incident field n* is the plane wave g*'^"^ '^ or the 
point source ^o{-,zo) with zq G CIq then / = — u^I^q, g = —du^ du\s(,, p = 0, q = 0, and if the 
incident field is the point source <I>i(-,zi) with zi G ^li then / = u^\so, 9 = duydu\sQ, p = 
-u'lsi, q = -du'/duls^. 

Theorem 2.3. The boundary value problem l\2.1\i — l\2.6^i admits at most one solution. 

Proof. Clearly, it is enough to show that u = in ^Iq, v = in 0,i and u; = in $72 for the 
corresponding homogeneous problem, that is, / = 5 = on 5o and p = q = on Si. Applying 
Green's first theorem over Qji, we obtain that 

du ^ f / ,2\ , f du 

u——ds= / (uAu+\\/u\ ) dx + u-—ds. 



dBn J So ^u 
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Using Green's first theorem over Jli and ^2 and taking into account the transmission conditions 
([231) and ([23]) . we have 

[ u^ds = I (uAn + |Vn|^) + Ao / [y /\v + \V v\^) dx 

JdBR Jur Jfli 

+A0A1 / {wAw + \Vw\'^)dx. (2.7) 

Making use of equations (j2.ip - (j2.3p and taking the imaginary part of ()2.7p we obtain, on noting 
that kQ, kf, /cl, Ao, Ai are positive numbers and Qn > 0, that 



/ u—ds = —k^XoXi'^ I n\w\^dx > 0. 
JdBR Jq, 



IdBR 

Thus, by Relhch's Lemma it follows that u = in W^\Bji. By the unique continuation 
principle, we have n = in Qq- Holmgren's uniqueness theorem [10] and the homogeneous 
transmission boundary conditions (j2.4p imply that v = in f^i. Finally, the transmission 
boundary conditions ()2.5p and the assumption on /cl give that w = in r22> which completes 
the proof of the theorem. □ 

For the proof of existence of solutions we need the fundamental solution of the Helmholtz 
equation with wave number kj [j = 0, 1, 2) given by 

= / x,2/GM3, x^y. (2.8) 

For z = 0, 1 and j = 0, 1, 2, define the single- and double-layer potentials Si^j and i^j.j, respec- 
tively, by 

(SijcP){x):= [ <^jix, y)<l){y)ds{y), x GR^S^, 

{K,^,<P){x) := ^ ^M^^^y)dsiy), x G R^\Si 

and the normal derivative operators K'- ■ and Tij by 



{K[^j4>){x) := -4"^ / ^j{x,y)(l){y)ds{y), x G R\Si, 

oi^{X) Jsi 



The restrictions on Si of these operators will be denoted by Sij, Kij, K[ ■ and Tjj- (z = 0, 1, j 
0, 1, 2), respectively. To prove the existence of solutions, we also need the volume potential 



{V4'){x) := ki / ^2{x,y)[n{y) - l](t){y)dy, x G 

and its normal derivative operator denoted by V . For mapping properties of these operators 
in the classical spaces of continuous and Holder continuous functions we refer the reader to the 
monographs of Colton and Kress [3^ 4j . 
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Theorem 2.4. The boundary value problem (jg.jp — l\2.6^ has a unique solution. Further, the 
solution satisfies the estimate 

) (2.9) 

for some positive constant C = C{a). 

Proof. Following [3j and [1] we seek the unique solution in the form 

u = Ao-fsTo.oV'o + 5'o,o(Ao in JIq, (2.10) 

V = i^ci'i/'o + 5*0,100 + AiKi,i-0i + 5'i,i(?;)i inl^i, (2.11) 
w = Ki^2'4^i + Si^2(t>i + Vw in 1^2 (2.12) 

with four densities G C^'"(So), (^o G C°'"(5o), G ^^'"(^i), € C°'"(5i). Then from the 
jump relations we see that the potentials n, v and w given by (j2.10p — ()2.12p solve the boundary 
value problem (12. ip — (12. 6p provided the densities satisfy the system of integral equations: 



tpQ + A(Ao-firo,o - -f^o,i)V'o + A(S'o,o - SQ^i)(j)Q - AAiKi^iV-i - AS'i,!^! = A/ on 5o, (2.13) 
+ AAo(ro,i - To,o)^o + A(Aoi^o,i - -^o,o)'/'o 

+AAoAiri,iV'i + AAoi?l,i</>i = -Ag on 5o, (2.14) 

V'l + fJ^Ko^iipo + fiSo,i(j)o 

+fj.{XiKi^i - Ki^2)4'i + fiiSi^i - Si^2)4>i - l^Vw = fip on Si, (2.15) 

+^Ai(ri,2 - Ti,i)^/;i + /i(Aiii'i 2 - K'i i)(t>i + fiXiV'w = -fiq on Si, (2.16) 

where A = 2/(Ao + 1) and ;U = 2/(Ai + 1). 
Define the product space 

X := C^'^'in^) X C7i'"(n7) X Ci'"(n^) X Ci'"(S'o) X C7°'"(So) x C^'"(5i) x C70'"(5i) 

and introduce the operator A : X X given in the matrix form: 
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^(Ai^l,2 - K,l 


y 



which is compact since all its entries are compact. Then the system (j2.10p — ()2.16p can be 
rewritten in the abbreviated form 

{I + A)U = R, (2.17) 

where U = {u, v, w, ipo, cfiQ, V'l, 4>i)'^ ^ R = (0) 0, 0, A/, —\g, fj.p, —fJ-q)^ and / is the identity opera- 
tor. Thus, the Riesz-Fredholm theory is applicable to establish the existence of solutions to the 
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system (j2.17p if we can prove the uniqueness of solutions to the system. To this end, let U be 
a solution of the homogeneous system corresponding to ()2.17p (that is, the system p.l7p with 
R = 0). Then it is enough to show that U = 0. 

We first prove that iPq = ^.q = on Sq. From the system (IZTTl) or (l2ln]l - (f2T6] ) with 
A/ = —\g = (since i? = 0) it is known that n, v, w defined in (j2.10p - (j2.12p satisfy the problem 
(12. ip — (12. 6p with f = g = 0. Thus, by the uniqueness Theorem 12.31 u = in f^o, v = m 
and w = in ^12- 

Now define 

V := Ko^iTpo + So,i(j)o + XiKi^iTpi + S'l,!^! in Qq, 

u := -Kofiipo - -^S'o.o'/'o in ^. 

Ao 

Then by the jump relations for single- and double-layer potentials we have 

y — v = iljo, — ti + u = on 5*0, (2.18) 
Ao 

dv dv du du 

^ = -(Po, ^ + Ao^ = -0o on 6o. (2.19) 

dv du dv dv 

Thus, V and u solve the homogeneous transmission problem 

Av + k^v = in Qq, An -|- kfu = in 

with the transmission conditions 

~ ~ ^ dv du 
V -u = 0, — = Xq— on bo- 
dv dv 

It is clear that v satisfies the radiation condition (j2.6p . By [31 Lemma 3.40], v = in ^Iq and 
u = in Consequently, from the boundary conditions ()2.18p and (j2.19p we conclude that 
^/'o = 00 = on ^Q. 

In a completely similar manner, we can also prove that ipi = 4>i = on Si. Thus, we 
have established the injectivity of the operator I + A and, by the Riesz-Fredholm theory, (/ -|- 
A)~'^ exists and is bounded in X. The estimate ()2.9p follows from (I2.17p . The proof is thus 
complete. □ 

We now establish the following mixed reciprocity relation which is needed for the inverse 
problem. 

Lemma 2.5. ( Mixed reciprocity relation.) For the scattering of plane waves u*(-, d) with d ^ S 
and point-sources ^{-^z) from the obstacle we have 



( 1 ^ 

—u^{z, — x), z G Qq, X G S, 

47r 

— v(z, —x), z € Oi, X & S. 

K 47r 



$°°(x, z) 



Proof. By Green's second theorem and the radiation condition (j2.6p we get 
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for z £ QqU Qi and x G 5. By [T, Theorem 2.5], we have the representation 

for z G fio U J^i and X G 5. Adding ([2:20]) to (p:2T]) gives 

for z G U r^i and x G S". 

We first consider the case z G CIq. Use the boundary condition on Sq and Green's second 
theorem to deduce that for z G f^o, x G 5, 

Now Green's second theorem gives that for z G and x G 5 
From Green's formula (see Theorem 2.4 in |4j) it follows that 

for z £ Qq and x G 5. Adding (j2.23p to the equation (j2.24p we deduce, with the help of the 
transmission condition on 5o and Green's second theorem, that for z G Qq and x £ S, 



By (I2.22P and ()2.25p together with the transmission condition on 5*1 and the equation ()2.3p we 
deduce that 
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AqAi 
AqAi 
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47r 
1 

47r 



dw{y,-x) dw{y,z) , , . . 

^i^' ^) — s"^^ — ^(y' o / X 



z)Aw{y, -x) - w{y, -x)Aw{y, z))ds{y) = 0. 

^2 



We now consider the case z E Oi. Using the transmission condition on So and Si, the 
equations (|2.2p and (j2.3p and Green's second theorem, we obtain that 
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Ao 
47r 
Ao 
A-K 
Ao 
47r 

A-K 



u^y,z)'-^^ll^-uiy,-x)^^]dsiy) 



du{y) 



du{y) 



11 X , ^ / dv{y,-x) d{v{y,z) +^i{y,z))\ 
[v{y, z) + $i(y, z)) ^(y, -x) — ds{y) 



dv{y) 



dv{y) 



(I X , ^ / ..dv{y,-x) d{v{y,z) +^i{y,z)y 
{v{y, z) + $i(y, z)) t>(y, -x) — ds{y) 



dv{y) 

^ / ^dv{y,-x) d^i{y,z)\ 



dv[y) 



dv{y) 



dv(%j) 



A, 



47r 



So 



7f MJ'ily.- 



dv{y) 
dv{y, —x) 



Air 



Si 



7f Mfi(y>- 



du{y) 
dv{y, -x) 



v{y, -x) 



dv{y) 
dv{y) 



ds{y) 



dv{y) 



Ao / 

where use has been made of Green's formula (see Theorem 2.1 in [4J 



□ 



Denote by D any of Or, f^i, ^2-, Sq, Si. Let x* G -D be an arbitrarily fixed point and let us 
introduce the space Cq{D) which consists of all continuous functions h G C{D\{x*}) with the 
property that 

lim \{x — x*)h[x)\ 

exists. It can be easily verify that Co(D) is a Banach space equipped with the weighted maximum 
norm 

ll^llo, d:= sup \{x — x*)h{x)\. 

x^x* , xGD 

To prove the unique determination of 5*0 in the inverse problem in the next section, we 
need to consider the the boundary value problem (|2.ip — (j2.6p with / = zo)\sq, g = 

—d^o{-,ZQ)/diy\so, p = 0, q = 0, where zq G Qq, and investigate the behavior of the solution v 
on some part of 5*0. 

Lemma 2.6. For given functions f G C^'"(S'o) and g G C^'^^Sq), assume that u G C^^Qq) fl 
C^^'^in^), V G C^i^i) n C^^'^ijh) and w G C72((72) n Ci'"(II^) are the solution of the problem 
(jg.jp — (jg.gp . Let X* G Sq and let Bi, B2 {B2 H = 0) be two small halls with center x* and 
radii ri, r2, respectively, satisfying that ri < r2- Then there exists a constant C > such that 

dv 

< C'dl/llo, So + ll^llo, So + ll/lll, a, So\Bi + \\g\\o, a, So\Bi) (2.26) 

oo, So\B2 



loo, So\B2 + 



du 
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Proof. We consider the system ()2.17p of integral equations derived from the boundary value 
problem ()2.ip — (12. 6p . In addition to the space X, we also consider the weighted spaces Cq := 
C7o(n^) X Co{Th) X Ci' "(H^) X C7o(5o) x C7o(5o) x Ci'"(5i) x The matrix operator ^ 

is also compact in Co since all entries of A are compact (see [HI E] ) . Prom the proof of Theorem 
2.41 we know that the operator I + A has a trivial null space in X. Therefore, by the Fredholm 
alternative applied to the dual system {X, Cq) with the bilinear form, the adjoint operator 
I + A' has a trivial null space in Cq. By the Fredholm alternative again, but now applied to the 
dual system (Cq, Cq) with the bilinear form, the operator I + A also has a trivial null space 
in Cq. Hence, by the Riesz- Fredholm theory, the system (|2.17p is also uniquely solvable in Cq, 
and the solution depends continuously on the right-hand side: 



\U\ 



< 



C(ll/l|o, 5o + Ibllo, So)- 
In particular, this implies that 

ll'^lloo, So\B2 ^ C'dl/llo, So + ll^llo, So)- 



olio, So + ll'Aollo, So + 11^1 111, a, Si + ll'/'lllo, a, Si 

(2.27) 



(2.28) 



Before proceeding to estimate dv/dv we establish the following estimate in the spaces of Holder 
continuous functions for (n, u;, -00) '/'Oi V^i; i?^!)- 

||f^||l, Q, 3 := ll'"llo, a, Qfl\_B3 + ll^llo, o, T7i\B3 + ll^lll, a, fl2 

+ IIV'o||l, a, So\B3 + ll</'o||o, a, So\Bs + IIV'llll, a, Si + ||</'l||o, a, Si 
< C'(||/||o,So + lbllo,So + ll/lll,a,So\Bi + ll5llo,a,So\Bi) (2-29) 

where is a ball of radius and centered at x* with ri < < r2. 

We now choose a function pi G C^(S'o) such that pi{x) = for x € Sq\B2 and pi{x) = 1 in 
the neighborhood of We also choose another function p2 G C'^{Sq) such that p2{x) = 1 for 
X G Sq\B3 and P2{x) = in the neighborhood of Bi. 

Split U up in the form 



U 



piv 
w 

Pi 00 

V 01 / 



+ 



({1-Pi)u\ 
{I- pi)v 
w 

(1 - /9i)^o 
(1 - pi)(j)Q 

Ipl 

01 



Up, + Ui^ 



pi 



and for a matrix W use Wp to denote the same matrix but with its first, second, fourth and 
fifth rows multiplied by p{x). Then it follows from ()2.17p that 



P2 



R 



P2 



Ap2 Up-, Ap2 Ui—p, . 



(2.30) 



The mapping operator U Ap^Up, is bounded from Cq into X since its kernel vanishes in a 
neighborhood of the diagonal x = y. Furthermore, we apply Theorems 2.30 and 2.31 in [3] and 
the analogous results for Tq^i, Kq ^,Ti^i,K-^ ^ to obtain that 

II^P2^1-P1 llo, a < C'll^f^l-pi llo, a < C||t^l-pi lloo < C||f^l|0) 
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where the norms ||C/||o, a and ||f/||oo are defined as follows: the first, second, fourth and fifth 
components of U are defined by the corresponding norms and its third, sixth and seventh 
components are equipped with the C^' "(02); C^' and C"' °'{Si) norms, respectively. From 

(fOOj) and (12:271) it is derived that 



\U\ 



0, a, 3 



< 
< 



ll, a, n2 

+ ll^o||o, a, So\B3 + \\4>o\\o, a, 5o\B3 + 11^1 IK, a, Si + 
C||t^P2llo, a < C{\\Rp2\\o,a + \\U\\o) 

C(ll/llo, a, So\Bi + 115110, a, So\Bi + ll/llo, So + Ibllo, So)- 



'lllO, a, Si 



(2.31) 



It remains to estimate ||^o||i, a, So\B3- Multiplying (j2.13p by P2{x) we obtain, on using (j2.30p 
and noting the fact that the integral operators mapping C'^'"-functions into C^'"-functions are 
bounded, that 

llV'olll.a, So\B3 ^ ll/'2V'o||l,«, So 

< C(||p2(Ao-f^O,0 - ^0,l)V'o||l,Q + \\P2{So,0 - 5'o,l)(/'o||l,a 
+ ||p2-?^l,lV'l||l,a + Il/O25'l,l0l||l,a + \\p2f\\l,a) 

< C(||C/||o + ||(l-pi)C/||o,a+||p2/||l,a) 

< C(||[/||o + ||C/||o,a,3 + ||/||l,a,So\Bj (2-32) 

where we have used the fact that ipQ = piTpQ + {1 — pi jipQ and (pQ = pi^o + (1 — Pi)</'o in order to 
get the third inequality. Combining ()2.27p and (|2.3ip - (|2.32p yields the desired estimate (|2.29p . 

Qy 

. From ()2.1ip and the jump relation we have that on So, 

0, a, So\B2 



We now estimate 



du 



dv 
dv 



1 



(2.33) 



Writing -00 = /CiV'o + (1 — /Oi)V'o and 0o = Pi4>q + (1 — Pi)4>o-, it follows from (|2.33p that 



dv 



du 



(1 -Pi TT 



< 

0, a, So\B2 "^^ 0, a, So 

< C(||C/||o + ||(l-/9l)V'o||l,a,So + 
+ a, Si + HlWo, a, Si) 

< C(||?7||o + ||t/||l,a,3). 



Pl)(P0\\0, a, So 



Combining this with ([2:271) and ([2:29]) gives 

dv 



dv 



0, a, So\B2 

This completes the proof. 



< C'(ll/llco(So) + llsllcoCSo) + ll/lll,a,Sa\Bi + II5'IIo,q,So\Bi ! 



(2.34) 



□ 



To prove the unique determination of Si, we consider the boundary value problem p.ip — ()2.6p 
in the case when the incident field is the point source ^i{-,zi) with zi G Oi, that is, / = 
''^'Iso) 9 = du^dvlso, p = —u^lsi, q = —dv}/dv\si- Arguing similarly as in the proof of Lemma 
12. 6[ we can obtain the following a priori estimate of the solution w on some part of 5*1. 
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Lemma 2.7. Given four functions f G C^'"(5q), 9 G C^'^iSo), P G C^'''{Si)_and q G C^'^'iSi), 
let u € C2(rJo) n Ci'"(no), G ^^(rJi) n ^^'"(ni) ond i« G C^fOa) n (^^'"(Qs) i/ie solution 
of the problem (fOj) - (fg7g|) . Lei x* G 5i and /ei Si, B2 (-B2 n = 0) 6e too smaZZ balls with 
center x* and radii ri, r2 (ri < r2), respectively. Then there exists a constant C > such that 



\w 



loo, Si\B2 + 



di 



< C'(||/||i,„, 5o + llffllo, a, So 

00, Si\B2 

+ \\P\\0, Si + Ikllo, Si + \\p\\l, a, Si\Bi + Ikllo, a, Si\i?J 



3 Unique determination of the interfaces 5*0 and 5*1 

Following the ideas of Kirsch and Kress [8j for transmission problems in a homogeneous medium, 
Kirsch and Paivarinta [9] for transmission problems in an inhomogeneous medium and Liu and 
Zhang [11] for impenetrable obstacle scattering in a piecewise homogeneous medium, we prove, in 
this section, that the interfaces Sq and Si can be uniquely determined by the far field pattern. To 

do this, we need the following four lemmas, in which $7 = ]R^\r2o so that C O and Q, = M^\r2o 
for some domain 0,q with interface Sq = O^Iq D dO, and domain 1^2 C 57. 



Lemma 3.1. For $72 C 57, 5^2 C 17 let Gq be the unbounded component ofW^ \ (17 U 57) and let 

d) = ■u°°(x, d) for all x, d G S, where u°°{x, d) is the far field pattern of the scattered field 
u'^{x,d) corresponding to the obstacle 572, the interface Sq and the same incident plane wave 
u^{x,d). For z G Gq let {u'^{-, z),v{-, z),w{-, z)) be the unique solution of the problem 

Au' + k^u'' = in57o\{z}, (3.1) 
Av + kjv = in57i, (3.2) 
Avu + k\nw = in 572, (3-3) 

- V = -<i>Q{-,z), AoTr = ^ — - on5o, 3.4 

v-w = 0, — = Ai— on5i, (3.5) 
ou 0: ■ 



V 



lim r { — - ikou' ] = 0. (3.6) 



Assume that {u^ {■, z),v{-, z),w{-, z)) is the unique solution of the problem (jg.jp — (jg.gjl with 
5^0, 57i, 572, 5*0, 5*1, n replaced by Q,q,^Ii,Q2, Sq, Si,n, respectively. Then we have 

u"^ {x, z) = V{x, z), X € Gq. 

Remark 3.2. By Theorem 12.41 the problem (j3.ip - (|3.6p has a unique solution. 
Proof. By Rellich's lemma [3], the assumption u°°{x, d) = u°°{x, d) for all x, d S implies that 

u'^(x, d) =u'^{x,d), x G Gq, d G S. 



Then, for the far field pattern corresponding to incident point-sources we have by Lemma fT. 
that 

$°°((i,z) = $°°(d,z), zeGQ,deS. 
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Thus, Rellich's lemma [1] implies that 

u^{x, z) = u^{x, z), X € Gq. 



□ 



Arguing similarly as in the proof of Lemma 13.11 we have the corresponding result for the 
point source in the bounded homogeneous medium. 



Lemma 3.3. For ^2 C ft, ft2 C ft let d := fl\ (J^2 Ufi2) and let u'^{x,d) = u°°{x,d) for 
all X, d G S where u°°{x,d) is the far field pattern of the scattered field u^{x,d) corresponding 
to the obstacle fl2, the interface Si and the same incident plane wave u^{x,d). For z G Gi let 
{u'^{-, z),v{-, z),w{-, z)) be the unique solution of the problem 







= 


in Oq, 


(3.7) 




Av + kfv = 


= 


in fli \ {z}, 


(3.8) 




Aw + k^nw = 


= 


in r22, 


(3.9) 


-V = ^i{x,z), — 

ov 


dv d^i{x, 

Ao a ~ a 
ov ov 




on 5o, 


(3.10) 


^ / ^ dv 

V - w = -9i{x,z), A 

ov 


dw d^i{x, 
^ dv dv 


z) 


on 5i, 


(3.11) 


lim r ( 

1 — ►oo \ 


du" \ 
-Q^ - ^^on^J = 


0. 




(3.12) 


V{-,z),v{-,z),w{-,z)) is the 


unique solution 


of the 


problem (3.1) — 


I\8.PA with 



fli,fl2, Si,n replaced by fti,ft2, Si,n, respectively. Then we have 

v{x, z) = v{x, z), X G Gi. 

Lemma 3.4. Assume that f £ L'^{fli),p G C{So),q G C{Si), g,r} e C{So) with ?/ / and 
T] < on Sq. Then the following problem has a unique solution v G C'^{fli) H C{fli) and 
w G C^{fl2)riG{Ih): 



Av + kfv = f 

+ k\nw = 
dv 

+ ir]v = g 

ov 



v — w 



p, 



dv ^ dw 
dv ^ dv 



Furthermore, there exists a constant C > such that 

Halloo. Til 



in fli, 
in ^2, 

on So, 
on Si. 

D, Si + Iklloo, 5i)- 



(3.13) 
(3.14) 

(3.15) 
(3.16) 



Proof. We first prove the uniqueness result, that is, f = in Qq, = in ili if / = in Qq, 
g = on Sq, p = q = on Si. With the help of equations (|3.13p - (j3.14p and boundary conditions 
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(I3l3]l - (l3l^ . we have 

{{Av + kfu)v} dx 



[ (-\Vv\'^ + k'^\v\'^)dx + [ v^ds- f v^ds 
Jn, Jso du Js, dv 

I {—\Vv\^ + k\\v\^)dx — i I ri\v'^ds — \i l w-^ds 
J J So J Si 



'Si du 

= / {-\Vv\'^ + kl\v\'^)dx - i r]\v\'^ds-Xi {-kln\w\'^ + \Vw\'^)dx. 

J J So ^ ^2 

Taking the imaginary part of the above equation gives that w = on some part F of Sq since 
both 7^ and r/ < on and Ai is a positive number. By the boundary condition (|3.15|) it 
follows that V = dvjdv = on F. Thus, by Holmgren's uniqueness theorem [lOj, v = in f^i. 
Using the transmission boundary conditions (|3.16|) and the assumption on k\^ we conclude that 
i(; = in r22. 

We now prove the existence of solutions using the integral equation method. To this end, 
introduce the volume potential 



{V*j){x) := [ '^i{x,y)f{y)dy, x e n^, 



which defines a bounded operator V* : L^(J7i) — > H'^(Qi) (see Theorem 8.2 in [3]). Now look 
for a solution in the form 

v{x)= -F*/ + ^o,i(0i + '/'2) + AiKi,i(V'i+V2) + ^i,i(xi + X2) in fii, (3.17) 
w{x)= T/*w; + i?i,2(^i + V'2) + ^i,2(xi + X2) in 1^2 (3.18) 

with six unknown densities 0i E C(S'o), </>2 G H^I^Sq), tpi G C(S'i), -02 £ H^{Si), xi S C(5i), 
X2 G H2(^Si). Then from the jump relations we see that the potentials v,w given by (j3.17p and 
(j3.18p solve the boundary value problem (|3.13p — (|3.16p provided the six densities satisfy the 
following system of integral equations: 

01 + 2(i^o,i + ivSo,i)<Pi + 2Ai(fi,i + iriKi^){^i + V2) + 2{k[^^ + ir?5i,i)(xi + X2) = 2g on 5o, 

<A2 + 2{k'^ + ir^So,i)(t>2 = 2 f |- + iri] {V* f) on So, 



du 

■01 - ijlV*w + /i5o,i(0i + 4)2) + p.{\iKi^i - Ki^2)tpi + ^(^1,1 - 5'i,2)xi = MP on Si 

ip2 + //(Aii^i^i - i^i,2)V'2 + - 5'i,2)x2 = l-iV*f on Si 

d ~ I II 

Xi + '"'^i^^*^ ~ f^^oA^i + 4'2) + AiAi(ri_2 - ri,i)V'i + ti{XiK^ 2 - -^i,i)xi = -m on Si 



I I d 

X2 + /i(AiA'^ 2 - ^l,l)X2 + MAi(ri,2 - ri,i)02 = ~i^Q^^^*f^ 



where /i = 2/(Ai + 1). This system is well defined since both the single and double layer 
potentials with //a -density are continuous and the restriction to the boundary of the derivative 
of the volume potential with L^-density is in H2. The trick to divide the density in a layer 
potential into a continuous and i/2-part has been used previously by Kirsch and Paivarinta [9]. 
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Precisely, we seek a solution {v,w,cj}i,(l)2,XijX2,Vi,V'2) £ Y := C{Qi) x C{Q2) x C{So) x 
H2[So) X C{Si) X H2(Si) X C(S'i) x H2(^Si) to the above system of eight integral equations. 
Standard arguments using the uniqueness of solutions to the problem show that this system 
admits at most one solution in Y. Therefore, the Riesz-Fredholm theorem is applicable and 
yields the existence of solutions to the boundary value problem with the estimate 



loo, Qi 



7=7 < 



c^(ll^*/lloo,n, + lblloo,so + 



d . 



{V*f) 



+ \\p\\oo,Si + Iklloo, Si) 
< + llflloo.So + Iblloo.Si + lkl|oo,5i; 

for some constant C > 0. The lemma is thus proved. 



□ 



In order to prove the unique determination of Si, we need the following lemma which was 
proved in [9] (see [9l Lemma 4.4]; note that our definition of the parameter rj and our assumption 
that is not a Neumann eigenvalue of Aw + k2nw = in are different from those in [9l 
Lemma 4.4]). 

Lemma 3.5. Assume that f € L'^{Q2) o-nd g,ri C{Si) with f] ^ and rj < on Sq. Then the 
following problem has a unique solution w € C'^{Qi) Ci C(J7i) of 



Aw + /c|n?i; 
dw 



f 



di 



+ irjw = g 



in 1^2, 
on Si. 



Furthermore, there exists a constant C > such that 

ll^lloo.na ^ <^(ll/llL2(ni) + IbllocSa)- 
We are now in a position to state and prove the main result of this section. 

Theorem 3.6. Assume that Aq ^ 1, \i ^ 1, let Sq, Sq be two penetrable interfaces and let 
Q2, ^2 be two penetrable obstacles for the corresponding scattering problem. If the far field 
patterns of the scattered fields for the same incident plane wave u*(x) = e^^^^'"^ coincide at a 
fixed frequency for all incident direction d S and observation direction x S, then Sq = Sq, 
Si = Si. 

Proof. We first prove that Sq = Sq. Let Gq be defined as in Lemma [3TT1 Assume that Sq ^ Sq. 
Then we may assume without loss of generality that there exists zq € Sq\Q. Let B2 be a small 
ball centered at zq such that B2 (1 Q = 9. Choose h > such that the sequence 



Zj :-- 



h , , 

Zq + -v{Zq), 

3 



1,2,..., 



is contained in Gq R B2, where v(^zq) is the outward normal to 5*0 at zq. Using the notations in 
Lemma [3?T] and letting (u?, w,, w^^ and ■Ui, i/;,) be the solutions of (I3.1I) - (I3.6D with z = z 

Then, by Lemma [3.H 



and {Uj, Vj, Wj) be the solutions of (j3.ip - (j3.6p with ^ — ^j. 
>j in Gq. Since zq has a positive distance from Q, we conclude 



from the well-posedness of the direct scattering problem that there exists C > such that 

duj 



\Uj\\oo, 5onB2 + 



du 



< G for all j > 1. 



(3.19) 



00, SonB2 
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Choose a small ball Bi with center zq which is strictly contained in B2. Since 



||^o(-,^i)||o,5o + ll^o(-,^i)||l,a,So\Bi < 
||^(-,2i)||o,So + ll^(-,^j)llo,a,So\Bi < C 

for some positive constant C independent of j, we conclude from Lemma 

dvj 



that 



IPjIloo, So\B2 + 
From this it follows that 



du 



< C for all j > 1. 



00, 5o\B2 



\XoVj - ^>o(-,2;j)||oo, So\B2 < C for ah 3 > 1, 



A, 



dvj d<^>o{-,Zj] 



< C 



00, 5o\B2 



du du 
The transmission boundary conditions yield 

\\vj - ^oi-,Zj)\\oo, SonB2 = \\uj\\oo, SonB2 < C 



for all i > 1. 



An 



du 







duj 




du 


00, SonB2 


du 


00, SoniJ2 



Combining (f3:2T]) and I^Jl^ yields 
Ao 



dvj d'^o{-,Zj] 



< C 



< C for all 3 > 1. 



for all J > 1, 
for all J > 1. 



00, So 



(3.20) 
(3.21) 

(3.22) 
(3.23) 

(3.24) 
(3.25) 



du du 

This can be used together with (j3.20p to prove the estimate 

\\XoVj - ^oi-,Zj)\\oo,So < C for all j > 1. 

In fact, choose a function r/ G C^(5o) supported in So\B2 with r] ^ and r/ < on S'q. Then 
V* := AoWj — $o('i -^j) solves the boundary value problem: 

Av* + klv* = {k^-kj)M;Zj) 
Aw* + k'^nw* = 

... . dw* _ dM;z,) 

du du du 

2 



in 




in 




on 


So 


on 


Si 



where w* = XoWj. Since, by (|CTD and (fX^ . / := (feg - A;f)$o(-,2i) G -^^^^i), 5 := 
(a/5i/ + ir?) [Ao^;, -^'o(-,^i)] G C(5o),p := -«'o(-,^i) G C{Si) andg := -a$o(-, ^i)/5i^ G C(5i), 
then the desired result ()3.25p follows from Lemma 13.41 

Now the triangle inequality together with (|3.22p and (j3.25p implies that 

||(Ao - l)^o(-,2:i)lloo,5onB2 < ||Ao^'o(-,^;j) - Aot'j||oo,SonB2 + IIAo^^j - ^o(-, 2;i)||oo,5onB2 

< C. 

This is a contradiction since Aq / 1 and l<I>o(2;o, ^j)|oo,Son-B2 ^ 00 as j — > 00. Thus, 5*0 = Sq. 

Arguing similarly as above, but using the point source $1 in Gi and Lemmas 12. 7^ 13.31 and 
13.51 we can easily prove that 5i = 5i. The proof is thus complete. □ 
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Remark 3.7. Our method can be extended straightforwardly to both the 2D case and the case 
of a multilayered medium, and a similar result can be obtained (that is, all the interfaces between 
the layered media as well as the shape of the embedded obstacle can be uniquely determined). 



4 Unique determination of the refractive index n 

In this section we shall prove a uniqueness theorem for recovering the refractive index n from 
the far field pattern. In doing this, we need the following completeness result. 

Lemma 4.1. The set {dw{-, d) /du\ d S} of normal derivatives of the fields w{-, d) correspond- 
ing to incident plane waves with directions d S are complete in L^(S'i) provided k"^ is not a 
Neumann eigenvalue of Aw + /clraw = in $72- 

Proof. Let Br be a large ball that contains $7 and such that fcg is not a Dirichlet eigenvalue of 
Alt; + A:|nw = in Br. Then the restriction to SBr of the set of plane waves e^^o^-d]^ d (z S is 
complete in L^{dBji) (see Theorem 5.5 in Thus we only have to show that the operator 
A : L^{dBji) — > L^(S'i) has a dense range Acj) = dw/du, where w solves the boundary problem 
()l.ip - (ll.6p and (j) is the boundary data of the interior Dirichlet problem 

An* + fcpn* = in Br, = 4> on OBr. 

A simple repeated application of Green's formulas yields that the L^-adjoint A* of A is given by 

A*V := 



f du* 


du \ 











where [u* ,v* ,w*) solves 





An* + klu* 


= 


in Qq, 


(4.1) 




Av* + klv* 


= 


in ill. 


(4.2) 




Aw* + A;|nw;* 


= 


in ^22, 


(4.3) 


v* = 0, 


du* dv* 

-T, -^0^— 

du dv 


= 


on 5o, 


(4.4) 


= -V^, 


dv* ^ dw* 
du ^ du 


= 


on 5*1, 


(4.5) 


lim r 1 




= 


r = \x\, 


(4.6) 



w 



and n is a solution of the interior Dirichlet problem 

Au + kQU = in Br, u = u* on dBR. 

We just need to show that A* is injective. Let A*^ = 0. Then we have du* /du = du/dv and 
n* = n on dBR. Define 



v*, in ^'\Br, 

u, in Br 
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Then v is an entire solution to the Helmholtz equation Av + k^v = in (see [7|) satisfying 
the radiation condition and therefore must vanish identically in R^. Hence, ti* = in R^\i?/j. 
By the unique continuation principle, u* = in JIq- Using the transmission conditions (j4.4p and 
Homogren's uniqueness theorem we conclude that = in f^i and therefore v* = dv* jdv = 
on 5i. By the transmission conditions ()4.5p . = on 5i. Thus we have t;;* = in 0,2 since 
k2 is not a Neumann eigenvalue of Aw + k2nw = Q 'mQ.2. Therefore, = by the transmission 
conditions (j4.5p . This completes the proof. □ 

Based on this completeness result, we are able to prove the following orthogonality relation. 

Lemma 4.2. If the far field patterns for the refractive indices n and n coincide, then for any 
solution w G C^(02) of Aw + /c|ni(; = in O'f^d any solution w G (7^(02) of Aw + /c|n{(; = 
in ^2 we have the relation: 

(n — n)w'wdx = 0. (4.7) 

Proof. We first prove (j4.7p for the special case w = w{-, d). By Rellich's lemma and Holmgren's 
uniqueness theorem, we have w{-,d) = w{-,d) and dw{-,d) / dv = dw{-,d)/dv on Si. From this 
and the equations satisfied by w{x,d),w{-,d) and w we deduce that 



/ 



(n — n)w{x,d)w{x)dx = I [Aw{x,d) + k2nw{x,d)]w{x)dx 



02 '^2 J^2 

1 



2 ; {^A[w{x,d) — w{x,d)] + k2n[w{x,d) — w{x^d)]^ w{x)dx 
= -T2 / [w{x,d) -w{x,d)]{Aw + klnw)dx = (4.8) 

where use has been made of Green's first theorem in getting the third equation. 

To complete the proof we need to show that any general w can be approximated in Lp'{Q.2) 
by functions w{-,d), d & S. Let us assume the opposite, that is, the set {w{-,d)\ d £ S} is not 
dense in L^(Q2) sense in W := {w £ (7^(172)1 Aw + A;|n7i; = in ^2}. Then by Hahn-Banach 
theorem there would be an / € L'^{Q2) such that 

f{x)w{x,d)dx = (4.9) 

O2 

for w{x, d) with all d £ S, but for some w S W, 

/ f{x)w{x)dx / 0. (4.10) 

7^2 

Let u G H'^{Q.2) be a solution to the interior Neumann problem 

Au + klnu = f in O2, 7^=0 on ^i. (4.11) 

dv 

From (14. 9p and (14. lip it follows on using Green's first theorem that 



0= / f{x)w{x,d)dx = / [Au{x) + k2nu{x)]w{x,d)dx = / u — — ^ — ds. 
Jq.2 JQ.2 Jsi 
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This, together with Lemma l4.ll imphes that u = on 5i. Thus, by Green's first theorem we 
conclude that 

Jq2 J^2 J^2 

which contradicts ()4.10p . The proof is thus complete. □ 

It is well known that the products ww of solutions to A^t; + A;2?^w = and Aw + k'^niv = for 
two different refractive indices n and n is complete in L^(ri2). Such a result was first established 
by Sylvester and Uhlmann [T7] and simplified by Hahner [5]. The reader is also referred to the 
argument of Theorem 6.2 in a recent review paper by Uhlmann [18]. Using this completeness 
result and Lemma 14.21 we have in fact proved that n = n. 

Theorem 4.3. Assume that the interfaces Sj {j = 0, 1) are known and fe^ is not a Neumann 
eigenvalue of Aw + k'^nw = in 0,2- Then the inhomogeneity n is uniquely determined by the 
far field pattern u°°{x, d) for all x, d & S. 

Remark 4.4. From Lemma [4. II and the argument in the proof of Lemma [4. 2 1 it follows that the 
Cauchy data of the solution for two different refractive index n and n coincide. This result can 
also be extended straightforwardly to the two dimensional case. Recently, Bukhgeim |2] proved 
that the refractive index n can also be uniquely determined from the set of Cauchy data in the 
2D case. Therefore, Theorem 14.31 also holds for the two dimensional case. 

Combining Theorem 13.61 with Theorem 14.31 we have the following uniqueness result. 

Theorem 4.5. Assume that Xq (7^ 1), Ai (7^ 1), k^, ki and k2 are given positive numbers 
and that is not a Neumann eigenvalue of Aw + /c^nii; = in 0,2- Then the interfaces 
Sj {j = 0, 1) and the inhomogeneity n are uniquely determined by the far field pattern n°°(x, d) 
for all X, d S. 
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